Intrinsic Photoconductivity of Ultracold Fermions in Optical Lattices 
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We report on the first experimental observation of a persistent alternating photocurrent in an 
ultracold gas of fermionic atoms in an optical lattice. The dynamics is induced and sustained by 
an external harmonic confinement. We find a counterintuitively momentum-dependent oscillation 
frequency for excited particles and a fast decay of holes which we attribute to spatial trapping. Life- 
time measurements reveal a significant enhancement of particle-hole recombination with increasing 
interactions. 
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Photoconductivity describes the change of a materials 
conductivity following an excitation with incident pho- 
tons. If the photon energy is resonant with a band tran- 
sition, electrons are excited from the valence band to the 
conduction band and an initial insulator becomes con- 
ducting PQ. Today, photoconductivity is widely used in 
technological applications such as semiconductor photo- 
diodes and photoresistors. It also provides a powerful 
probe for novel materials, such as graphene [2], transis- 
tors made from carbon nanotubes [5] or semiconductor 
nanowires [3] . To extend the understanding of such com- 
plex materials, atomic quantum gases have proven to be 
powerful model systems. In this context, however, it is 
desirable to develop and adopt versatile probing meth- 
ods [5H7]. Owing to its excitational structure in several 
bands, photoconductivity can provide deeper insight into 
intra- and interband dynamics as well as orbital effects, 
which gained much interest in recent years. In the field 
of quantum gases, multiband interactions and dynamics 
have been experimentally studied mainly with bosonic 
atoms [5Hl7j. For fermionic atoms, population trans- 
fers into higher bands have been recently demonstrated 
via Landau-Zener transitions across Dirac cones [TS] and 
momentum-resolved lattice amplitude modulation [1.9;. 

In this letter, we study experimentally and theoret- 
ically the photoconductivity of fermionic atoms in an 
optical lattice. We create uncoupled particle and hole 
excitations using lattice amplitude modulation and thor- 
oughly investigate the dynamics of the atoms in the sec- 
ond excited band as well as the holes in the valence band. 
The lattice modulation corresponds to the incident pho- 
tons in solid state photoconductivity experiments, while 
the atoms (holes) correspond to the electrons (holes) in 
the conduction (valence) band. A typical difference be- 
tween ultracold gas experiments and real materials is the 
presence of an overall harmonic confinement, which - in 
our experiment - plays the role of an external force. It in- 
duces pronounced oscillations of the atoms in the excited 
band with a counterintuitive dispersion relation which 
differs maximally from the zero lattice case for shallow 
lattices. We explain this with a reduced effective mass in 
the conduction band. In strong contrast to the long-lived 
coherent dynamics of the atoms, we observe a dramati- 



cally shorter lifetime of the holes in the valence band. We 
attribute this effect - generic to harmonically trapped lat- 
tice systems - to spatial trapping at localized states at the 
edge of the atomic cloud. As a key result of this paper, 
photoconduction of fermionic atoms is mainly connected 
to long-lived populations in excited bands, forming the 
majority component, whereas holes in the valence band 
decay quickly. 

For the experiments we prepare either a spin-polarized 
Fermi gas with m = 9/2 or an interacting binary spin 
mixture of m — —9/2 and m = —5/2 in the / = 9/2 
groundstate manifold of 40 K (for details see [20). The 
atoms are held in the combined potential of an optical 
dipole trap with variable trapping frequency loq and an 
optical lattice at A = 1030 nm. To induce a photocurrent, 
we excite the system via lattice modulation spectroscopy 
[2TJ [22] as explained in [19]. This two-photon process 
creates particle excitations in the conduction band and 
leaves vacancies (holes) in the valence band. Both have 
the same initial quasimomentum qo- Due to the dif- 
ferent curvature of the involved bands, qo can be arbi- 
trarily tuned by choosing a particular modulation fre- 
quency [see Fig. [lja)]. To record the time evolution of 
the excited system, we detect the quasimomentum dis- 
tribution by performing adiabatic bandmapping followed 
by resonant absorption-imaging after 15 ms time-of- flight 
[T9] [23]. Recall, that the bandmapping technique maps 
particles in the second excited band to real momenta in 
the interval k/kBz — [2,3] and k/k BZ = [—3,-2], con- 
stituting the third Brioullin zone. Particles in the lowest 
band are mapped to fc/fcez = [— 1, 1], with fcez = 2tt/X. 

A typical time evolution for a spin-polarized gas is 
shown in Fig. [ljb) after an excitation at a quasimomen- 
tum qo = 0.5 and a lattice depth of s = 10 in units 
of the recoil energy E T = /i 2 fcg Z /2m with m the mass 
of K. The atoms in the conduction band show a pro- 
nounced oscillation in momentum space. The decay time 
of the excitations is of the order of 100 ms. This indicates 
a very slow recombination between particles and holes. 
However, the holes in the valence band undergo a fast de- 
cay and are only visible in momentum space for the first 
2 ms. This indicates a local trapping of the holes within 
the valence band, reminiscent of trapping of charge car- 
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FIG. 1. (a) Principle of excitation via lattice amplitude mod- 
ulation (dotted lines) and band mapping process (solid lines) 
in momentum space, (b) Measured photocurrent after excita- 
tion to the conduction band at 10-B r - Shown are the column 
densities of the momentum distribution at different times af- 
ter the excitation. Atoms in the valence band are represented 
by the central plateau. The excitations in the conduction 
band clearly oscillate in momentum space, (c) Center-of- 
mass quasimomentum of the particle-excitation in the con- 
duction band at 10i? r for two different trapping frequencies 
cjo = 2-7T x 66 Hz (circles) and cjrj = 2ir x 50 Hz (diamonds). 
Solid lines are fits to the data. Extracted oscillation frequen- 
cies are Q = 2n x 295 ± 9 Hz (circles) and ft = 2ir x 213 ± 3 Hz 
(diamonds) . 



riers in photoconducting solids [T]. Consequently, the 
excited atoms form the majority component of the pho- 
tocurrent in our system. We concentrate on their behav- 
ior in the following and afterwards analyze the trapping 
of the holes in more detail. Our method is related to the 
measurements in [13] using bosonic atoms, where oscilla- 
tions in the spatial domain were studied. Hence, neither 
holes nor excitations with small go could be investigated. 

We first examine the effect of the harmonic trapping 
potential on the dynamics of the particles. Figure [ljc) 
shows two typical measurements with different harmonic 
trapping frequencies ujq for elsewise equal initial param- 
eters. Fits to the data reveal, that the oscillation fre- 
quency linearly depends on uiq but is substantially higher 
than the bare trapping frequency. This result signifi- 
cantly deviates from the dynamics in harmonic poten- 
tials without any lattice, where the oscillation frequency 
is equal to ujq. To understand the main features of the 
photocurrent we first present two simplified descriptions 
before we describe the complete experimental and nu- 
merical results. 

Consider the case of a particle at a given quasimomen- 
tum in an excited band of an optical lattice. This is 
depicted in Fig. [2ja) using the extended zone scheme. In 
the presence of an harmonic confinement different quasi- 
momenta are coupled which causes an oscillation in mo- 
mentum space. The existence of a substantial band gap 
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FIG. 2. (a) Sketch of the time evolution of atoms in the con- 
duction band using the extended zone scheme, (b) Oscillation 
frequency in the conduction band at go = for different lat- 
tice depths, calculated with (1). Inset shows the complete 
conduction band for 4, 8, 12 and 16 E r . Note the the decrease 
in curvature with increasing lattice depth, leading to an in- 
crease of the effective mass m' . 



between first and second excited band leads to a striking 
modification of the simple harmonic oscillator dynamics: 
When a particle with initial quasimomentum q reaches 
the gap, it is Bragg reflected and thus omits a huge part 
of the extended zone scheme. The oscillation continues 
until the particle reaches —go and is reversed. Thus, de- 
pending on go, the intraband dynamics is limited to dif- 
ferent fractions of the harmonic oscillation. This leads 
to a strong momentum dependence of the oscillation fre- 
quency w(go). In summary, by omitting the inner part of 
the spectrum, the oscillation period is dramatically de- 
creased in comparison to the pure harmonic case. Recall 
that the oscillation is still driven by the harmonic trap 
and thus the frequency is proportional to ujq as observed 
experiment ally. 

For oscillations with small initial quasimomentum, 
there is a complementary explanation which exploits the 
form of the bands in the reduced zone scheme. It shows, 
that the photocurrent frequency non-trivially depends 
on the lattice depth: The presence of an optical lattice 
results in a modified kinetic energy with a momentum- 
dependent effective mass. Around q — 0, where the con- 
duction band has a minimum, this effective kinetic en- 
ergy can be approximated as a quadratic dispersion with 
a renormalized mass m! . The effective Hamiltonian is 
thus given by 
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where m is the real mass of the particles and w = 
^m/m'ujQ the renormalized trapping frequency. Espe- 
cially for shallow lattices the effective mass m! becomes 
much smaller than m and thus uj becomes very large, as 
shown in Fig.^b). This result leads to the counterintu- 
itive feature, that the dynamics in shallow lattices max- 
imally differs from the zero lattice case. This apparent 
contradiction is resolved by the onset of Landau-Zener 
tunneling to the first excited band in very shallow lat- 
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FIG. 3. (a) Comparison of measured frequencies and nu- 
merical results for different lattice depth and quasi-momenta. 
Filled symbols represent measurements with non-interacting 
mixtures. Open symbols show data with interacting mixtures. 
Lines show numerical calculations using the full trapping and 
lattice potential. Note the strong increase of w for shallow lat- 
tices and small q. (b) Exemplary comparison of experimental 
data (left) and numerical result (rigth) for the time evolution 
in the conduction band. Only the third Brioullin zone is de- 
picted, (c) Decay time of particles from the conduction band 
as a function of the scattering length for a 3D lattice of 8 i5 r . 
Shown is the time after the excitation until only 1/e atoms are 
still in the conduction band. All error bars solely correspond 
to fit errors, representing two standard deviations. 



tices. The bandgap becomes small and vanishes in the 
zero-lattice limit, where the atoms oscillate in the har- 
monic confinement with the bare trapping frequency u>q. 
In contrast to the conduction band, in the valence band 
the effective mass ml is always larger than m for any lat- 
tice depth and thus the oscillation frequency is smaller 
than ujq. Therefore, photoconductivity in the conduction 
band qualitatively differs from the conduction in the va- 
lence band. 

We measured the photocurrent in our system to in- 
vestigate all abovementioned dependencies. Figure [3^ a) 
shows our measurements of the photocurrent frequency 
as a function of the lattice depth and the initial quasi- 
momentum. For spin-polarized gases, we measured the 
frequency of the conduction band photocurrent at s = 8 
and s = 10 for different go- In the interacting case we 
investigated the three lattice depths s = 4, 8, 16 at the 
background scattering length of 168.5 ao [23]. All data 
exhibit the behavior expected from the argumentation 
presented in connection to Fig. [2j The oscillation fre- 
quency to decreases with increasing lattice depth and in- 
creasing q . 

To obtain a quantitative prediction for the dispersion 
relation u>(qo,s) of the photocurrent, we performed nu- 



merical calculations. We assume a single particle in the 
conduction band, represented by a Gaussian distribution 
of Bloch states centered around a given quasimomentum 
go- This state is subject to the full potential of the opti- 
cal lattice and the harmonic confinement. For further de- 
tails see [2D] . The results of our calculation are presented 
in Fig. [3|a) in comparison with our experimental data, 
showing very good agreement. Figure (3^b) displays an 
exemplary comparison of experimental data in momen- 
tum space and a corresponding calculation. The oscil- 
latory behavior with a renormalized frequency is clearly 
reproduced by the numerical results. In particular, the 
calculations confirm the experimentally determined lin- 
ear dependence between u> and the bare trapping fre- 
quency ui for all g . Moreover, we find that for small g , 
where ([!]) is valid, it agrees with the numerical results. 

For the interacting case, no significant deviations from 
the single-particle picture could be observed in the oscil- 
lation frequency. In contrast to the spin-polarized case, 
however, we find a substantially enhanced decay rate of 
atoms from the conduction band. This includes both, loss 
of particles from the trap and decay of particles into lower 
bands. To investigate the decay-dependence on the inter- 
action strength, we used a Feshbach resonance at 224 G 
[25] with a width of 7.6 G [26]. The corresponding data is 
shown in Fig. [3^c) . The total decay rate from the conduc- 
tion band strongly depends on the interaction, while the 
total atom loss is independent of the scattering length, 
which we checked independently. Thus, the results are 
a measure for the time needed for recombination of free 
charge carriers, which increases with decreasing interac- 
tion as naively expected and shows a maximum at van- 
ishing interaction. As an application, this effect might be 
used as a new probe to characterize more precisely the 
zero-crossing of Feshbach resonances. In conclusion, we 
find very good agreement of experiment and numerical 
calculations for many different parameters, which shows, 
that photoconductivity of excited fermions in optical lat- 
tice with superimposed harmonic traps can be thoroughly 
understood. 

In addition to the excited particles in the conduction 
band, one can clearly observe holes in the valence band 
in the photoconduction measurement of Fig. |TJb) . As 
shown in more detail in Fig. [4^a) , the holes decay how- 
ever, very fast on the timescale of a few ms, which cannot 
be explained by recombination with excited atoms, which 
have a much longer lifetime [see Fig. [jjjc)]. The decay 
of the holes in momentum space instead indicates spa- 
tial trapping, which is also observed in photoconducting 
normal insulators pQ. In these systems, local imperfec- 
tions of the periodic potential typically lead to trapping 
of the holes, which can no longer participate in the pho- 
tocurrent anymore. In ultracold atom experiments, it is 
well known, that the harmonic potential leads to local- 
ized states at the edge of the system [27]. These states 
are off-resonant for tunneling in the lattice, since the lo- 
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FIG. 4. (a) Time evolution of relative hole depth for different 
configurations. Circles: 4iJ r , go = 0.0, Wq = 20 Hz, interact- 
ing mixture; Diamonds: 4E r , qo = 0.0, ujo = 36 Hz, interact- 
ing mixture; Triangles: 1015 r , go = 0.5, luo = 50 Hz; Squares: 
6E T , qo — 0.4, cjo = 63 Hz. Solid lines are single-particle cal- 
culations as described in the text, (b) Exemplary comparison 
of experimental data (left) and numerical result (right) for the 
time evolution in the valence band. Only the first Brioullin 
zone is depicted, (c) Relative hole depth after 250 fj,s for dif- 
ferent trapping frequencies at 4iJ r , go = 0.5 in an interacting 
mixture. Solid line is a single-particle calculation with an oc- 
cupation of 40 states in the lattice direction. All error bars 
solely correspond to fit errors, representing two standard de- 
viations, (d) Numerical results for the relative hole depth as 
in (c) for uio = 2tt x 50 Hz (dashed) and 2n x 90 Hz (solid) 
as a function of the number of occupied states in the lattice 
direction. 

cal potential difference between two neighboring sites is 
larger than the bandwidth of the lattice. This is true 
both for particles and for holes. However, for a hole to 
be dynamically trapped in localized states, in analogy to 
real solids, it is necessary, that these state are occupied 
by particles before the excitation. Since the localized 
states are at the edge of the atomic cloud, their occupa- 
tion increases with the filling and the lifetime of the holes 
should decrease accordingly. 

To quantitatively analyze the trapping of the holes in 
our experiment, we adopt a theoretical description from 
solid state physics, where holes in an otherwise com- 
pletely filled valence band can be described as parti- 
cles with negative mass. We have simulated the holes 
in analogy to the excited particles by assuming a single 
particle with negative mass in the valence band of the 
combined periodic and harmonic potential (see |20j for 
details). Figure |4^a) shows, that the decay of the holes 
can be very accurately described by our model for various 
regimes. Examplarily, in Fig.[4][b) a comparison of exper- 
imental and numerical data in momentum space shows 



the fast decay of a hole in less than 1 ms. 

As for the excited atoms in the conduction band, we 
investigated the influence of the harmonic trapping in 
detail. Figure [4^c) shows the decay of a hole in depen- 
dence of u>o- The lifetime of the holes clearly decreases 
with increasing ujq in good agreement with our numerical 
calculations. The main effect of the increased confine- 
ment is to increase the local potential difference between 
two neighboring sites. This creates more localized states 
near the center of the system. At constant total filling, 
the result is a higher relative population of such localized 
states. Assuming local trapping of the holes, this explains 
the decrease in lifetime. We also numerically checked the 
dependence on the filling at constant trapping strength. 
The results are shown in Fig. Qd) for two different trap- 
ping frequencies loq. At high filling the holes decay much 
faster than for a system with low filling. This supports 
our experimental findings from Fig. [4j^c) and our intuitive 
explanation. 

In conclusion we have presented a comprehensive study 
of photoconductivity in an ultracold gas of fermionic 
atoms. By independently analyzing the dynamics of ex- 
cited particles and holes we have shown that atoms in the 
conduction band constitute the majority charge carriers 
in our system. The observed long-lived oscillatory dy- 
namics could be reproduced very well by our numerical 
simulations and proves to be very sensitive to the har- 
monic as well as the lattice potential. In particular we 
measure that counterintuitively the dynamics maximally 
differs from the purely harmonically trapped case in very 
shallow lattices. Opposed to the persistent dynamics 
of the atoms in the excited band, we observe a very 
short free-carrier lifetime of the holes which we attribute 
to trapping in localized states at the edge of the com- 
bined lattice and harmonic potential. These results may 
prove crucial for further studies on particle-hole excita- 
tions such as excitons. The presented measurements ex- 
tend the available techniques to explore dynamical prop- 
erties of optical lattice systems and equally important 
emphasize the increased role of the harmonic confinement 
for experiments involvong excited spatial bands. Finally, 
our results constitute an important contribution for the 
understanding of fundamental dynamical properties of 
fermionic quantum gases in optical lattices. 
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SUPPLEMENTAL INFORMATION 

This supplemental information discusses the 
preparation of our atomic sample (SI), details of 
the fitting procedures for the experimental data 
(S2,S3) and the theoretical calculations both for 
particle (S4) and hole (S5) excitations. 

SI. PREPARATION OF THE ATOMIC SAMPLE 

By sympathetic cooling, we create a mixture of spin- 
polarized 87 Rb and 40 K atoms in a magnetic trap. The 
atoms are transferred adiabatically to a crossed optical 
dipole trap operated at 811 nm with a 1/e 2 radius of 
120 /im. After switching off the magnetic trap, we re- 
move the rubidium atoms from the trap using a reso- 
nant light pulse. For experiments with spin polarized 
atoms, the preparation is finished her. For experiments 
with interacting mixtures, we use a series of rf-pulses 
and -sweeps to prepare an equal mixture of the hyperfine 
states m = —9/2 and m = —5/2 in the hyperfine mani- 
fold / = 9/2. This mixture is evaporatively cooled in the 
following by reducing the laser power of the optical dipole 
trap. The final particle number is about N = 5-10 4 atoms 
at typical temperatures of 0.2 Tp. After the preparation 
we linearly ramp up an optical lattice within 100 ms. The 
lattice consists of up to three orthogonal retro-reflected 
laser beams at A = 1030 nm with a 1/e 2 radius of 200 /xm. 
For measurements at the Feshbach resonance, the mag- 
netic field was set to the final value 50 ms prior the 100 ms 
optical lattice ramp. 

To initialize the photocurrent, we modulate the am- 
plitude of one of the lattice directions for 1 ms with a 
frequency, that is resonant with a transition from the 
lowest energy band to the second excited band. This ex- 
cites a fraction of particles and leaves vacancies in the 
lowest band. Due to the different curvature of the bands, 
the resonance condition depends on the quasimomentum. 
By tuning the modulation frequency, we have full control 
over the quasimomentum of the excited particles. Since 
lattice amplitude modulation does not imprint any quasi- 
momentum, the holes in the lowest energy band have the 
same quasimomentum as the particles. 

S2. ANALYSIS OF EXPERIMENTAL DATA: 
PARTICLE EXCITATIONS 

In this section we describe how we extract the oscilla- 
tion frequency of the excited atoms from the experimen- 
tal data. For each time step we determine the quasimo- 
mentum of the excitation by taking the center-of-masses 
of the atoms in the conduction band independently at 
positive and negative momentum. To be insensitive to 
global displacements, we measure the difference of both 



excitation centers instead of their absolute positions. We 
extract the oscillation frequencies f2 from the differen- 
tial center-of-masses by fitting an exponentially damped 
cosine of the form 

Aq(t) = A exp(-rt) cos(fit + $) + C , 

with oscillation amplitude A, damping rate T, a phase 
shift <f> and a constant offset C. In the experiment, al- 
ways two excitations with opposite initial quasimomcnta 
<7o and — q^ are created. Both excitations independently 
perform oscillations in the combined potential of lattice 
and harmonic trap. After a quarter of an oscillation 
period, both excitations arrive at q = 0, are Bragg re- 
flected and continue at the other side of the Brioullin 
zone, respectively. Since the excitations are not distin- 
guishable, our center of mass determination cannot re- 
solve this. Therefore the extracted data artificially ex- 
hibits a turning point of the oscillations at this position. 
The data thus shows an oscillation with twice the funda- 
mental frequency: = 2u>. Consequently, the factor of 
2 is corrected in all experimental data except for Fig. 1, 
such that we always show to instead of fl 

S3. ANALYSIS OF EXPERIMENTAL DATA: 
HOLE EXCITATIONS 

To determine the depth of the hole excitation in the 
valence band, we fit a sum of three Gaussians to the first 
Brioullin zone in the time-of-fiight picture. One of the 
Gaussians represents the atomic background density in 
the valence band. The other two represent the holes at 
positive and negative q, respectively. The background is 
determined from the momentum distribution directly af- 
ter the excitation for each time series individually. For 
the dynamical evolution we take the form of the back- 
ground to be constant for all times and solely allow for a 
variation of the absolute magnitude as a function of time. 

S4. CALCULATION OF CONDUCTION BAND 
DYNAMICS 

In this section we describe the numerical calculations 
for the particle excitations in the conduction band. To 
derive quantitative predictions for the oscillations for all 
initial quasi-momenta, we performed calculations includ- 
ing both the periodic and the harmonic confinement. We 
assume a single particle confined to the full potential 

p 2 1 
H = 2m + sEr cos ^ BzX "> 2 + 2 mw ° x2 ' 

By diagonalizing the homogeneous lattice Hamiltonian 
i?o = p 2 /2m + sE r cos{k bzx) 2 for a single particle, we 
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obtain the Bloch states for a given lattice depth s. In 
energy space, the system has bands of allowed states di- 
vided by gaps where no states are located. In the ex- 
periment, the initial excitation is produced by a lattice 
amplitude modulation pulse of t = 1 ms duration. This 
pulse width is much larger than the trapping frequency 
ujo- For the parameters discussed here, it holds in gen- 
eral, that 

E c » h/t » Huj , 

where Ec is the band width of the conduction band. 
Therefore, the harmonic confinement is negligible dur- 
ing the preparation. In all calculations we assume a 
distribution of Bloch states centered at a certain quasi- 
momentum q . For simplicity we use a Gaussian dis- 
tribution, with a variance a = 364.5 Hz, corresponding 
to the width of the 1 ms lattice modulation. We obtain 
the time evolution by exact diagonalization of the full 
Hamiltonian H which leads to a non-trivial dynamics as 
shown in Fig. 3(b). To extract the oscillation frequen- 
cies from the dynamics, we calculate the center-of-mass 
of the quasi-momentum distribution with respect to the 
homogeneous lattice Hamiltonian H and fit a cosine to 
the data. For a typical calculation we use up to 400 
quasi-momcnta and 11 bands. 

In the calculations we only use one excitation which 
has a positive quasi-momentum. We clearly observe the 
Bragg reflection at q = and obtain an oscillation fre- 
quency in agreement with the approximation of equa- 
tion (1). As mentioned in section S2, the analysis of 



the experimental data results in twice the oscillation fre- 
quency, since excitations at positive and negative quasi- 
momentum are indistinguishable. This effect is consis- 
tently reproduced by the calculations, if we incoherently 
overlap two excitations with opposite sign in the quasi- 
momentum, as shown in Fig. 3(b). 



S5. CALCULATION OF HOLE DYNAMICS 

For electron gases in solids a missing electron can be 
described as a single particle with negative mass. Such a 
hole is in complete analogy to an electron. We adopted 
this prescription for holes in the valence band of our har- 
monically trapped quantum gas. To describe this situa- 
tion numerically, we assumed a single particle with neg- 
ative mass m* — ~m in the valence band and calculated 
its evolution in the presence of the harmonic confine- 
ment and the lattice potential. The excitation has the 
same shape as in the conduction band, since the missing 
atoms in the valence band directly correspond to excited 
atoms in the conduction band. We also take into account 
the finite filling of the valence band. This is done by us- 
ing only the lowest energy states in the time evolution of 
the initial state. Finally, we include finite temperature 
by using a Fermi-Dirac distribution instead of a sharp 
edge when weighting the included eigenstates of H. To 
match the occupation number to the experimental situ- 
ation, we take the total atom number and the trapping 
frequencies and calculate the number of occupied states 
in the appropriate direction. 



